Chapter 13 

SOC Systems in Astrophysics 



by Markus J. Aschwanden 

The universe is full of nonlinear energy dissipation processes, which occur intermit- 
tently, triggered by local instabilities, and can be understood in terms of the self- 
organized criticality (SOC) concept. In Table 2.1 we included a number of cosmic pro- 
cesses with SOC behavior On the largest scale, galaxy formation may be triggered by 
gravitational collapses (at least in the top-down scenario), which form concentrations of 
stars in spiral-like structures due the conservation of the angular momentum. Similarly, 
stars and planets form randomly by local gravitational collapses of interstellar molecu- 
lar clouds. Blazars (blazing quasi-stellar objects) are active galactic nuclei that have a 
special geometry with their relativistic jets pointed towards the Earth, producing erratic 
bursts of synchrotron radiation in radio and X-rays. Soft gamma repeaters are strongly 
magnetized neutron stars that produce crust quakes (in analogy to earthquakes) caused 
by magnetic stresses and star crust fractures. Similarly, pulsars emit giant pulses of 
radio and hard X-ray bursts during time glitches of their otherwise very periodic pulsar 
signal. Blackhole objects are believed to emit erratic pulses by magnetic instabilities 
created in the accretion disk due to rotational shear motion. Cosmic ray particles are 
the result of a long-lasting series of particle acceleration processes accumulated inside 
and outside of our galaxy, which is manifested in a powerlaw-like energy spectrum 
extending over more than 10 orders of magnitude. Solar and stellar flares are produced 
by magnetic reconnection processes, which are observed as impulsive bursts in many 
wavelengths. Also phenomena in our solar system exhibit powerlaw-like size distri- 
butions, such as Saturn ring particles, asteroids, or lunar craters, which are believed 
to be generated by collisional fragmentation processes or their consequences (in form 
of meteroid impacts). The magnetosphere of planets spawns magnetic reconnection 
processes also, giving rise to substorms and auroras. 

While these astrophysical processes have been interpreted in terms of the self- 
organized criticality concept (for a comprehensive overview see Aschwanden 201 la), 
quantitative theoretical modeling of astrophysical SOC phenomena is still largely un- 
explored. In the following Section 13.1 we outline a general theory approach to SOC 
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phenomena, which consists of a universal (physics-free) mathematical/statistical as- 
pect, as well as a physical aspect that is unique to each astrophysical SOC phenomenon 
or observed wavelength range. In the subsequent Section 13.2 we discuss then the as- 
trophysical observations and compare the observed size distributions with the theoreti- 
cally predicted ones. 

13.1 Theory 

A system with nonUnear energy dissipation governed by self-organized criticaUty (SOC) 
is usually modeled by means of cellular automaton (CA) simulations (BTW model; 
Bak et al. 1987, 1988). A theoretical definition of a SOC system thus can be derived 
from the mathematical rules of a CA algorithm, which includes: (1) A S-dimensional 
rectangular lattice grid, (2) a place-holder for a physical quantity Zij^k associated with 
each cellular node Xij^k, (3) a definition of a critical threshold Zcrii, (4) a random input 
Azij^k in space and time; (5) a mathematical re-distribution rule that is appUed when 
a local physical quantity exceeds the critical threshold value which adjusts the state of 
the next-neighbor cells, and (6) iterative time steps to update the system state Zij^k{t) 
as a function of time t. Although this definition is sufficient to set up a numerical 
simulation that mimics the dynamical behavior of a SOC system, it does not quantify 
the resulting powerlaw-like size distributions in an explicit way, nor does is include 
any physical scahng law that is involved in the relationship between statistical SOC 
parameters and astrophysical observables. A quantitative SOC theory should be gen- 
eralized in such a way that it encompasses both the mathematical/statistical aspects of 
a SOC system, as well as the physical scaling laws between observables and statistical 
SOC parameters. In the following we generaUze the fractal-diffusive SOC model (FD- 
SOC), described in Aschwanden (2012a) and outlined in Section 2.2.2, which includes 
three essential parts: two universal statistical aspects, i.e., (i) the scale-free probability 
theorem, (ii) the fractal-diffusive spatio-temporal relationship, and a physical aspect, 
i.e., (iii) physical scaling laws between geometric SOC parameters and astrophysical 
observables, which may be different for each observed SOC phenomenon and each ob- 
served wavelength in astrophysical data. Some basic examples of physical scahng laws 
are derived for fragmentation processes, for thermal emission of astrophysical plasmas, 
and for astrophysical particle acceleration mechanisms. 

13.1.1 The Scale-Free Probability Theorem 

Powerlaw-like size distributions are an omnipresent manifestation of SOC phenomena, 
a property that is also called a "scale-free" parameter distribution, because no pre- 
ferred scale is singled out by the process. Of course, the scale-free parameter range, 
over which a size distribution exhibits a powerlaw function, is always limited by in- 
strumental sensitivity or a detection threshold at the lower end, and by the finite length 
of the time duration over which a SOC system is observed and sampled, at the upper 
end. Bak et al. (1987, 1988) associated the scale-free behavior with the fundamental 
property of 1/f-noise that is omnipresent in many physical systems, giving rise to a 
power spectrum of P{i') oc 
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Figure 13.1: The geometric relationships between the Eudidean avalanche volvmie Ve oc L^, 
the fractal avalanche volume V (x L^' , and the system volume Vsys oc L^ys is visualized in 3D 
space {S = 3), as a function of the avalanche length scale L and system size Lsy^. Note that the 
probability for an avalanche with size L scales with the ratio of the system volume Vsys to the 
EucUdean avalanche volume Ve. 



However, here we give a more elementary explanation for the powerlaw behavior 
of SOC size distributions, namely in terms of the statistical probability for scale-free 
avalanche size distributions. A key property of SOC avalanches is that random distur- 
bances can produce both small-scale as well as unpredictable large-scale avalanches of 
any size, within the hmitations of a finite system size Lgys at the upper end, and some 
"atomic" graininess AL at the lower end (i.e., a sand grain in sand avalanches, or the 
spatial pixel size AL in a computer lattice grid). The statistical probabihty disttibution 
N{L) for avalanches with size L can be calculated from the statistical probability. If 
no particular size is preferred in a scale-free process, the number N{L) of possible 
avalanches in a S-dimensional system with a volume Vgys = Lf^^ is simply the system 
volume Vsys divided by the Euclidean volume Ve = of a single avalanche with 
length scale L (Fig. 13.1), 

N{L) oc (^^^ oc L-^ , (13.1) 

In a slowly-driven SOC system, none or only one avalanche happens at one particu- 
lar time, but the relative probability that an avalanche can happen still scales with the 
reciprocal volume. This is our first basic assumption of our generalized SOC model, 
which we call the scale-free probability theorem. This theorem itself predicts a power- 
law function for the basic size distribution of spatial scales, and is distinctly different 
from the gaussian distribution function that results from binomial probabilities. For 
3-dimensional SOC phenomena (S = 3), thus we expect a size distribution or differ- 
ential occurrence frequency distribution of N{L) oc L~^, or a cumulative occurrence 
frequency distribution of A^(> L) oc L^^. 

We can now derive the expected size distributions for related geometric parameters, 
such as for the Euclidean avalanche area Ae or volume Ve. If we simply define the 
Euclidean area Ae in terms of the squared length scale, i.e., Ae oc or L oc A^J"^, 
which has the derivative dL/dAE oc A^ ' , we obtain for the area size distribution 
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N{Ae), using N{L) oc L"^ (Eq. 13.1), 



N{AE)dAE oc N[L{Ae) 



dL 



dAi 



dAE oc A-^^+^^/^rf^E , 



(13.2) 



yielding N{Ae) oc A^ for 3D phenomena (S = 3). Similarly we define the Eu- 
cUdean volume, i.e., Ve oc L"^ or L oc V^^^, which has the derivative dL/dVE oc 
V^/^~'^, yielding a volume size distribution NiVE), 



N{VE)dVE oc N[L{Ve)] 



dL 



dVE 



dVE oc V'^^-^'^UVe , 



(13.3) 



yielding NiVE) oc Vg for 3D phenomena (S" = 3). 

In the case that avalanche volumes are fractal, such as characterized with a Haus- 
dorff dimension Ds in EucUdean space with dimension S, the fractal volume V scales 
as, 

y oc , (13.4) 

which yields L oc V^/^^ and the derivative dL/dV oc ^jj^ tjjus the size 

distribution. 



7V(y)dy oc N[L{V)] 



dL 



dV 



(13.5) 



The Euclidean limit of non-fractal avalanches would yield for 5 = 3 and Ds = S = 3 
the same powerlaw index ay = 1 + {S — 1)/Ds = 5/3 as derived for in Eq. (13.3). 
For fractal avalanches, with Dg « (1 + S) /2 = 2.0 for 5 = 3, we obtain a slightly 
steeper powerlaw distribution, N{V) oc than for Euchdean avalanches, i.e., 

N{Ve) cx Ve'^^ 

A similar effect occurs for fractal avalanche areas A. If we assume an fractal struc- 
ture with Hausdorff dimension D2 in 2D space 



(13.6) 



which yields L oc A^^^^ and the derivative dL/dA = A^^^^ ^, and thus a size 
distribution of 



N{A)dA oc N[L{A)] 



dL 



dA 



dA oc 



(13.7) 



For fractal avalanches, with D2 ~ (1 + S')/2 = 1.5 for the 2-D fractal dimension, 
we obtain a powerlaw distribution N{A) oc which is sUghtly steeper than for 

Euclidean avalanche areas, i.e., N{Ae) oc A~^^'^. 

In some astrophysical observations, such as in solar and planetary physics, the areas 
A of SOC phenomena can be measured, while spatially integrated emission (or spa- 
tially unresolved emission) in (optically-thin) soft X-ray or extreme-ultraviolet wave- 
lengths is often roughly proportional to the fractal volume V of the emitting source, 
and thus the size distributions Eqs. (13.5) and (13.7) can be used to test our scale-free 
probability theorem (Eq. 13.1). 
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If the scale-free probability theorem (Eq. 13.1) is correct, the derived size distribu- 
tions for spatial scales N{L), avalanche areas N{A), and avalanche volumes N{V), 
should be universally valid for SOC phenomena, without any physical scaUng laws. 
They should be equally valid for earthquakes or solar flares, regardless of the physi- 
cal mechanism that is involved in the nonhnear energy dissipation process of a SOC 
event. In Section 13.2 we will present some astrophysical measurements of size distri- 
butions of such geometric parameters (L, A) which can corroborate our assumption of 
the scale-free probability theorem. Volume parameters (V) can usually not directly be 
measured for astrophysical objects, except by means of stereoscopy or tomography of 
nearby objects. 

13.1.2 The Fractal-Diffusive Spatio-Temporal Relationship 

After we have established a framework for the statistics of spatial or geometric pa- 
rameters of SOC avalanche events, we turn now to temporal parameters, which can be 
defined by a spatio-temporal relationship. The temporal evolution of SOC avalanches is 
governed by the complexity of next-neighbor interactions above some threshold value, 
which has an erratically fluctuating time characteristics according to cellular automaton 
simulations. However, the mean radius r{t) of an evolving SOC avalanche was found 
to closely mimic a time dependence of r{t) oc t^^^, which can be associated with a 
classical random-walk or diffusion process (Aschwanden 2012a). Measurements of 
the spatial evolution of solar flares, which are considered to be an established SOC 
phenomenon, revealed a similar evolution, but tend to be sub-diffusive for the analyzed 
dataset (Aschwanden 2012b). Moreover, the instantaneous avalanche area was found to 
have a fractal structure, while the time-integrated avalanche area is nearly space-filling 
and thus can be described with an Euclidean area or volume. These two properties of 
fractal geometry and diffusive evolution have been combined in the fractal-diffusive 
SOC avalanche model (FD-SOC) (Aschwanden 2012a). We generaUze this concept 
now also for anomalous diffusion, 

r(t) oc K , (13.8) 

where k is the diffusion coefficient and the diffusive index /3 combines classical diffu- 
sion (/? = 1), as well as anomlous diffusion (/3 7^ 1). Anomalous diffusion processes 
include both sub-diffusion (/3 = G...1), as well as super-diffusion (^ = 1...2), also 
called hyper-diffusion or Levy flights, 

{/3 < 1 (sub-diffusion) 
P = 1 (classical diffusion) (13.9) 
/3 > 1 (super-diffusion or Levy flights) 

We show the generic time evolution of a sub-diffusion process with /3 = 1/2 and a 
super-diffusion process with /? = 3/2 in Fig. 13.2. Anomalous diffusion implies more 
complex properties of the diffusive medium than a homogeneous structure, which may 
include an inhomogeneous fluid or fractal properties of the diffusive medium. 

The spatio-temporal evolution of an instabiUty generally starts with an exponential 
growth phase (which we may call the acceleration phase), followed by a saturation or 
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Figure 13.2: Comparison of spatio-temporal evolution models: Logistic growth with parame- 
ters th = 1-0, Too = 1.0, TG ~ 0.1, sub-diffusion (/3 = 1/2), classical diffusion {j3 = 1), Levy 
flights or super-diffusion (/? = 3/2), and linear expansion (r oc t). All curves intersect att = tL 
and have the same speed v = (dr/dt) at the intersection point at time t = tL- (Aschwanden 
2012b). 



quenching phase (which we may call deceleration phase). In the logistic growth model 
(Section 2.2.1), the deceleration phase saturates asymptotically at a fixed value, while 
diffusive models do not converge but slow down progressively with time (Fig. 13.2). 

The diffusive scaling (Eq. 13.9) impUes then also a statistical correlation between 
spatial L and temporal scales T, 

L oc , (13.10) 

where T is the time duration of a SOC avalanche. From the scale-free probabihty the- 
orem N{L) oc (Eq. 13.1) we can then directly compute the expected occurrence 
frequency distribution N{T) for time durations T, 



N{T)dT oc N[L{T)] 



dL 



dT 



dT oc T-[i+(^-iW21 dT . (13.11) 



For instance, for 3D SOC phenomena {S = 3) we expect a powerlaw distribution 
N{T) oc r-(i+''), which amounts to N{T) oc T'^ o for classical diffusion, N{T) oc 
2^-1.5 Jqj. ^ sub-diffusion case {(3 — 0.5), or N(T) oc T~^-^ for a super-diffusion case 
(/3 — 1.5). The exponentially growing phase is neglected in this derivation, which 
implies a slight underestimate of the number of short time scales. Time scales can also 
directly be measured in most SOC phenomena, and thus provide an immediat test of 
the fractal-diffusive assumption made here, regardless of the physical process that is 
involved in the observed signal of SOC avalanches. 
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13.1.3 Size Distributions of Astrophysical Observables 

The previous theory on geometric (L, A. V) and temporal (T) parameters should be 
universally valid for the statistics of SOC phenomena, and thus constitutes a purely 
"physics-free" mathematical or statistical property of SOC systems. All other observ- 
ables of SOC events, however, are related to a physical (nonlinear) energy dissipation 
process, which needs to be modeled in terms of a correlation or scaling law with respect 
to the physics-free spatio-temporal SOC parameters. Say, if we observe a physical SOC 
variable x that has a correlation or powerlaw scaling law of x oc with the geometric 
SOC parameter L, we can infer the expected size distribution N{x)dx by substituting 
the scaling law. 

What is most common in astrophysical observations is the flux F or intensity that 
is observed in some wavlength range A (with physical units of energy per time), origi- 
nating from a source with unknown volume V. The flux F{t) can exhibit strong fluc- 
tuations during an energy dissipation event, but we can characterize the time profile of 
the event with a peak flux P, or with the time-integrated flux, also called fluence (with 
physical untis of energy), which we may denote with E. For optically-thin emission 
observed in soft X-ray or EUV emission, the emissivity or flux is approximately pro- 
portional to the source volume V, so it is most useful to quantify a scaling law of the 
flux F with the 3D volume V, which we characterize with a powerlaw index 7, 

FocV^. (13.12) 

From the size distribution of the fractal volume N{V) oc T/-[i+('5-i)/Os] (Eq. 13.5) 
and the scaling law V oi F'^/"' (Eq. 13.12) and its derivative dV/dF oc i^'^/f ~^ we can 
then derive the size distribution N{F) of fluxes F, 



N{F)dF oc N[V{F)] 



dV 



dF 



dF oc . (13.13) 



which has a typical powerlaw index of ai? ~ 2.0 (for S = 3, Ds « {1 + S)/2 2.0, 
and 7 f» 1). For peak fluxes P we have the same distribution, except for the fractal 
dimension having its maximum Euclidean value Ds ~ S, which yields. 



N{P)dP oc N[V{P)] 



dV 



dP 



dP oc p-[i+('5-i)/7S]^j^ _ 



which has a typical value of ap « 5/3 « 1.67. Finally, the total flux or fluence 
E = J F{t)dtK, FT, is found to have a size distribution of. 



N{E)dE oc N[V{E)] 



dV 



dE 



dE oc £-[i+(5-i)/(70s+2//3)]^^ _ (13,15) 



which has a typical value of ue ~ 3/2 = 1.5, for S = 3, Ds ^ {I + S)/2 k. 2.0, 
7 « 1, and /? « 1. 

In summary, if we denote the occurrence frequency distributions N(x) of a param- 
eter X with a powerlaw distribution with power index a^, 

N{x)dx oc a;-"- dx , (13.16) 
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OiL 


= S 




aA 


= 1 + {S- 


1)/D2 


av 


= 1 + {S~ 




aT 


= 1 + {S- 


l)P/2 


ap 


= 1 + {S- 


l)/(7^s) 


ap 


= l + (5- 


l)/(7^) 


ap 


= 1 + (5 - 


l)/(7£>s + 2//3) 



we have the following powerlaw coefficients for the parameters x = L,A,V,T,F,P, 
and E, 



(13.17) 



Thus, the various powerlaw indices depend on four fundamental parameters: the Eu- 
chdean dimension S of the SOC system, the fractal dimension Dg of SOC avalanches, 
the diffusion index /3 of the SOC avalanche evolution, and the scaUng law index 7 
between the observed flux and the SOC avalanche volume. Note, that the powerlaw 
slopes of the geometric (L, A, V) and flux parameters {F, P) do not depend on the dif- 
fusion index P, and thus are identical for classical or anomalous diffusion. Only the 
powerlaw slopes of the length scale, time scale, and peak flux (L, T, P) do not depend 
on the fractal dimension. All flux-related parameters (F, P, E) depend on a physical 
scaling law (7), which may be different for every observed wavelength range. 

In the following we generally assume 3D SOC phenomena (S = 3), for which the 
fractal dimension can be estimated by the mean value between the minimum and max- 
imum dimension where SOC avalanches can propagate coherently via next-neighbor 
interactions (which limits the minimum fractal dimension to Ds,min ~ 1-0 and maxi- 
mum fractal dimension to Dmax = S), 

,min ~\~ fjiax 

» 2 = — 2 — ' (13.18) 

which yields k. 2.0 for S = 3 and simplifies the powerlaw indices to 

aL =2 

aA = 7/3 

av =2 

ap =l + P^2 . (13.19) 

ap = 1 + 1/7 « 2 

ap = 1 + 2/(37) ~ 5/3 

ap = 1 + 1/(7+ l/^)« 3/2 

In astrophysics, the distributions of geometric parameters (L, A, V) can only be de- 
termined from imaging observations with sufficient spatial resolution (in magneto- 
spheric, hehospheric, and solar physics), while the distributions of aU other parameters 
(T, F, P, and E) can be measured from any non-imaging observations, such as from 
point-like stellar objects. 



13.1.4 Scaling Laws for Thermal Emission of Astrophysical Plas- 
mas 

Solar flares and stellar flares are observed in soft X-ray and extreme-ultraviolet (EUV) 
wavelengths, where the observed intensity is measured in a particular wavelength range 
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A given by the instrumental filter response function. Soft X-ray and EUV emission is 
produced by photons via the free-free bremsstrahlung process, free-bound transitions, 
or radiative recombination. In strong magnetic fields, cyclotron and gyrosynchrotron 
emission is also produced at radio wavelengths. Soft X-ray and EUV emission occur 
usually in the optically thin regime, and thus the total emission measure EM, which 
is proportional to the observed intensity in a given wavelength A, is proportional to the 
volume V of the emitting source. 

Fx oc EM = j n2(x) dV oc (n^) V . (13.20) 

Thus, if the electron density in a source would be constant, or the same among 
different flare events, we would have just the simple relation Fx oc (Eq. 13.12) with 
the scaling law index 7 = 1, which is an approximation that is often made. In fact, this 
is quite a reasonable approximation for measurements with a narrow-band temperature 
filter, which is sensitive to a particular electron temperature Tg, and thus probes also a 
particular range of electron densities and plasma pressure p = 2nekBTe that depend 
on this temperature Tg, whatever the scaUng law between electron temperature Tg and 
electron density n,. is. 

The proportionality constant between the flux intensity Fx and emission measure 
EM is dependent on the wavelength range A, because each wavelength filter is cen- 
tered around a different temperature range Tg that corresponds to the line formation 
temperature in the observed wavelengths A. Physical scaling laws have been derived 
to quantify the relationship between electron temperature Tg, electron density rig, and 
the spatial length scale Ljoop of coronal loops, e.g., by assuming a balance between the 
heating rate, conductive, and radiative loss rate (i.e., the so-called RTV law; Rosner, 
Tucker, and Vaiana 1978), being (in cgs-units), 

Te » 1400 {j)Lioopfl^ . (13.21) 

which we can express in terms of the electron density rig, using the definition of the 
ideal gas law, p = 2nekBTe, 

Thus, if there is no particular correlation between the loop length Lioop of the densest 
flare loops (with the highest emission measure) and the volume V of the active region, 
the density is only a function of the electron temperature, oc Tg . Consequently, 
if a narrowband temperature filter is used in a soft X-ray or EUV wavelength range, 
sensitive to a peak temperature Tx, the corresponding electron density is given in a 
narrow range also, rig oc T|, and thus the flux is essentially proportional to the flare 
volume (according to Eq. 13.20), Fx oc V, with a scaling index of 7 « 1. On the 
other hand, if a different powerlaw index 7 ^ 1 is measured, such an observation would 
reveal a systematic scaling of some parameters (Tg, rig, Lioop) of the densest flare loops 
with the size L of the active region. 

Another important quantity we want to calculate is the size distribution of thermal 
energies Eth, for which we expect a scaling law of (using Eq. 13.22), 

Eth = SnekeTeV oc T^V/Li^op , (13.23) 



10 



CHAPTER 13. SOC SYSTEMS IN ASTROPHYSICS 



where the most dominant value of the electron temperature is given by the peak 
of the differential emission measure distribution (DEM). Observationally, it was found 
that the DEM peak temperature of a flare scales approximately with the size i of a 
flare, i.e., Te x L (Aschwanden 1999), which yields with V oc L^^ (assuming that 
Lioop with the highest emission measure is uncorrelated with the active region size L), 

Eth oc L^+^s ^ (13.24) 

or Eth oc for Dg ~ 2.0 and 5 = 3. The size distribution for thermal energies is 
then expected to be, 



N{Eth)dEth = N{L[Eth\) 



dL 



dEth 



dEth=E;^'+^'-'y^'+''^^^dEth (13.25) 



which yields N{Eth.) oc E'^^-'^ for 5 = 3 and Dg w 2.0. However, we have to 
be aware that this estimate of the thermal energy contained in a flare requires multi- 
thermal measurements to derive the peak DEM temperature and cannot be obtained 
from a single narrowband filter measurement. Note that the size distribution of thermal 
energies with powerlaw index ast^ = 1.4 is very similar to the total energy in photons 
in any wavelength range, i.e., fluence, as = 1-5 (Eq. 13.19). 

The foregoing model for the thermal energy requires a scaling law between the 
flare size L and its statistical temperature T^. In practice, however, thermal energies 
were often estimated in a limited temperature range from the filter ratio of two narrow- 
band filters. Such filter ratio measurements are sensitive to a particular temperature Tg 
and electron density Ue (Eq. 13.22), which are then essentially constants in the expres- 
sion for the thermal energy, and thus the thermal energy is mainly proportional to the 
volume, 

Eth,v oc V , (13.26) 

so that the size distribution of thermal energies V is identical to the size distribution of 
volumes V, which has a powerlaw slope of ath.v = oiv = 1 + (5 — 1) /-D3, yielding 
values in the range of ath,v = 1-67 — 2.0, depending on fractal (D^ « (1 + 5)/2 = 2 
for 5 = 3) or Euclidean (D3 = S = 3) volume measurements. 

In some studies, the volume is approximated with a "pill-box" geometry, i.e., the 
product of the measured flare area A with a constant height h along the line-of-sight, 
V = Ah, which makes the volume proportional to the area, and consequently the 
thermal energy is mainly proportional to the flare area A, 

Eth,A oc A , (13.27) 

leading to a size distribution of Eth,A that is identical with that of the area A, which 
has the powerlaw slope ath,A = = 1 + (S* — I)/ D2, yielding values in the range 
of a.th.A = 2.0 — 2.3, depending on fractal {D2 ~ (1 + S')/2 = 1.5) or Euchdean 
(D2 — 2) area measurements in 5 = 3 space. Thus, accurate measurements of these 
size distributions provide a sensitive tool to diagnose the underlying physical scaUng 
laws and model assumptions. 



13.2. OBSERVATIONS 



11 



13.1.5 Scaling Laws for Astrophysical Acceleration Mechanisms 

Let us consider some simple examples of particle acceleration processes in astrophys- 
ical plasmas, such as solar flares, stellar flares, or cosmic rays. The simplest parti- 
cle acceleration process is a coherent direct current (DC) electric field, which can be 
characterized by an acceleration constant a over a system length L. Newtonian (non- 
relativistic) mechanics predicts for an electron with mass nie a velocity v = at after 
a distance L = (l/2)at^, which corresponds to a kinetic energy El of (where the 
subscript L refers to the length scale of the accelerator), 

El = ^ruev'^ = rUeaL , (13.28) 

which implies a linear energy increase with system length. El oc L. Thus, the size 
distribution of energies El is identical with that of length scales L, which we obtain 
from the scale-free probability theorem N{L) oc (Eq. 13.1), 

dL 



N{EL)dEL = N{L[El]) 



dEL 



dEL = El^dEL , (13.29) 



yielding an energy spectrum N{El) ~ EJ^^ for a 3D Euclidean volume. This simplest 
case is given by coherent acceleration thoughout the entire source volume. 

Another approach could be made by assuming that the magnetic energy density 
dEs /dV = B^/Stt per volume element converted into nonthermal particle energy is 
uncorrelated with the system size L, in which case the total converted magnetic energy 
just scales with the Euclidean volume V (x of the energy dissipation process, 
Eb oc Ve oc L^. The corresponding size distribution of magnetic energy Eb is then 
expected to scale as. 



N{EB)dEB = N{L[Eb) 



dL 



dEi 

5/3, 



dEB = E~^'^~'^'^UEb , (13.30) 



which translates into iV(£'B) oc E^ for a 3D Euclidean volume (S=3). Interestingly, 
this value is also identical with the size distribution of the peak flux P. Alternative 
scaling laws using the Alfvenic crossing time through the flare volume have also been 
considered (e.g, Shibata and Yokoyama 1999, 2002; Nishizuka et al. 2008). 

Thus, the two acceleration models yield powerlaw size distributions in the range 
of Ueb = 1-67 to = 3.0. The measurement of size distributions of nonthermal 
energies thus can yield valuable diagnostics about the physical nature of the underlying 
particle acceleration process. 

A summary of all theoreticaUy derived powerlaw indices expected in astrophysical 
systems is compiled in Table 1. 



13.2 Observations 

We discuss now a number of astrophysical observations with regard to their observed 
size distributions, which we compare with the foregoing theoretical predictions. A 
more detailed review of these measurements is given in Aschwanden (201 1 ; chapters 7 
and 8). 
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Table 13.1: Summary of powerlaw indices predicted in astrophysical systems, as a function of 
the dimensionality S, the fractal dimension Ds, the diffusion power index /3, and the energy- 
volume scaling index 7. 



Parameter 


Powerlaw index 
(general expression) 


Powerlaw index for 

S = 3, Da = 2, D2 = 3/2 
/3= 1,7 = 1 


Length scale L 


aL = S 


az, = 3 


Area A 


aA = l + {S- l)/D2 


aA = 7/3 


Volume V 


av = l + {S- 1)1 D3 


av = 2 


Time duration T 


aT = l + h 


ar = 2 


Flux F 


a_F = 1 + 1/7 


ap = 2 


Peak flux P 


ap = 1 + 2/(37) 


ap = 5/3 


Fluence E 


as = 1 + 1/(7+1/^) 


ap = 3/2 


Emission measure EMx 


"EMa = ay 


apMx = 2 


Thermal energy 


aE,^ = l + {S-l)l{i + Ds) 


c^EtH = 7/5 


Thermal energy Eth,A 


OiEth,A = aA 


"Bth =7/3 


Thermal energy Eth,v 


asth.v = av 


aEtH = 2 


Linear energy El 


= S 


= 3 


Magnetic energy Ejj 


<^En = {'2- 1/5) 





13.2.1 Lunar Craters 

The size of lunar craters was measured from pictures recorded with the lunar orbiters 
Ranger 7, 8, 9 by Cross (1966). A size distribution of 1,600 lunar craters, sampled in 
the Mare Tranquillitatis using data from Ranger 8, within a range of 0.56 to 69,000 
m, is shown in Fig. 13.3, exhibiting a powerlaw distribution ranging over 5 orders of 
magnitude with a slope of ^ « 2.0 for the cumulative occurrence frequency distribu- 
tion, which translates into a powerlaw slope of a w /3 -|- 1 w 3.0 for the differential 
occurrence frequency distribution, 

N{L) oc L-^ . (13.31) 

This corresponds exactly to our prediction of the scale-free probability theorem for 
avalanche events in 3D-space (S = 3). A similar powerlaw index of ap = 2.75 
was also found for the size distribution of meteorites and space debris from man-made 
rockets and satelUtes (Fig. 3.1 1 in Sornette 2004). 

How can we interpret this result? The justification of our scale-free probability 
theorem is the fact that the relative probability of partitioning a system into smaller 
parts scales reciprocally with the volume (Fig. 13.1). The leading theory of lunar crater 
formation is that their origin was caused by impacts of meteorites, and thus the size 
distribution reflects that of the impacting meteors and meteorites, which probably were 
produced by numerous random collisions, similar to the origin of planets, asteroids, and 
planetesimals. Both the Moon and the Earth were subjected to intense bombardment 
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Figure 13.3: Left; The lunar crater Daedalus, about 93 km in diameter, was photographed by 
the crew of Apollo 11 as they orbited the Moon in 1969 (NASA photo AS 11-44-6611). Right: 
Cumulative frequency distribution of crater diameters measured from Ranger 8 in the lunar Mare 
Tranquillitatis (Cross 1966). 

of solar system bodies between 4.6 and 4.0 billion years ago, which was the final stage 
of the sweep-up of debris left over from the formation of the solar system. 

An alternative explanation of lunar or terrestrial craters is a volcanic origin. If 
we attribute volcanoes with nonlinear energy dissipation avalanches in a slowly-driven 
system of stressing planet crust motions and build-up of subtectonic lava pressure, 
volcanic eruptions can also be understood as SOC phenomena. 

13.2.2 Asteroid Belt 

The asteroid size distribution has been studied in the Palomar Leiden Survey (Van 
Houten et al. 1970) and Spacewatcti Surveys (Jedicke and Metcalfe 1998), where a 
power law of N'^^"^{> L) cx L~^-^ was found for the cumulative size distribution 
of larger asteroids (L > 5 km), which corresponds to a differential powerlaw slope of 
az, ~ 2.8. In a Sloan Digital Sky Survey collaboration (Fig. 13.4, right), a broken pow- 
erlaw was found with N{L) oc L~^-^ for large asteroids (5-50 km) and N{L) oc L~'^ 
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Figure 13.4: Left top: A picture of tiie near-Eartli asteroid Eros with a size of 30 km, pic- 
tured by a space probe. Left bottom: The main asteroid belt located between the Jupiter and 
Mars orbit. The subgroup of Trojan asteroids are leading and trailing along the Jupiter orbit. 
(Courtesy of NASA/Johns Hopkins University Applied Physics Laboratory). Right: Differential 
size distribution of asteroids observed in the Sloan Digital Sky Survey collaboration (Ivezic et 
al. 2001). 



for smaller asteroids (0.5-5 km) (Ivezic et al. 2001). In the Subaru Main-Belt Asteroid 
Survey, a cumulative size distribution Af™™(> L) oc 2^-i-29±o.o2 found for small 
asteroids with L « 0.6 — 1.0 km (Yoshida et al. 2003; Yoshida and Nakamura 2007), 
which corresponds to a differential powerlaw slope of sa 2.3. Thus, the observed 
range « 2.3 — 4.0 of the powerlaw slopes of length scales L is centered around 
the theoretically expected value of = 3.0, predicted by our scale-free probability 
theorem. 

The origin of asteroids is thought to be a left-over distribution of planetesimals 
during the formation of planets, which were either too small to form bigger planets by 
self-gravitation, or orbited in an unstable region of the solar system where Mars and 
Jupiter constantly cause gravitational disturbances that prevented the formation of an- 
other planet. Thus, the final distribution of the asteroid belt is likely to be influenced by 
both the primordial distribution of the solar system as well as by recent collisions and 
further fragmentation of planetesimals. The collisional fragmentation process can be 
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considered as a mechanical instability that occurs in a multi-body gravitational field. 
The coUisional process is self-organizing in the sense that the N-body celestial mechan- 
ics keeps the structure of the asteroid belt more or less stable, despite of the combined 
effects of self-gravity, gravitational disturbances, collisions, depletions, and captures 
of incoming new bodies. The quasi- stability of the asteroid belt warrants the critical 
threshold in form of a finite colUsion probability maintained by the proximity of the 
co-orbiting asteroid bodies. 

13.2.3 Saturn Ring 

The distribution of particle sizes in Saturn's ring was determined with radio occupation 
observations using data from the Voyager 1 spacecraft and a scattering model, which 
exhibited a powerlaw distribution of N{r) oc (Fig. 13.5) in the range of 1 mm 
< r < 20 m (Zebker et al. 1985; French and Nicholson 2000). The size distribution 
revealed sUghtly different powerlaw slopes in each ring zone, e.g., = 2.74 — 3.03 
for ring A, = 2 . 79 for the Cassini division, or = 3.05 — 3.22 for ring C (Zebker 
et al. 1985). These results, again, are consistent with a fragmentation process that obeys 
the scale-free probabiUty theorem, similar to the distribution of sizes of asteroids and 
lunar craters, and predicts a size distribution of N{L) oc L~^. The conclusion that 
Saturn's ring particles are formed from the (collisional) breakup of larger particles, 
rather than from original condensation as small particles, was already raised earlier 
(Greenberget al. 1977). 

The Saturn ring, which consists of particles ranging in size from 1 cm to 10 m, is 
located at a distance of 7,000-80,000 km above Saturn's equator, and has a total mass 
of 3 X 10^^ kg, just about a little less than the moon Mimas. The origin of the ring 
is believed to come either from leftover material of the formation of Saturn itself, or 
from the tidal disruptions of a former moon. The celestial mechanics of the Saturn 
rings is quite complex, revealing numerous gaps in orbits that have a harmonic ratio 
in their period with one of the 62 (confirmed) moons (with 13 moons having a size 
larger than 50 km). Similar to the asteroid belt, the Saturn ring can be considered as 
a self-organizing system in the sense that the gravity of Saturn and the gravitational 
disturbances caused by Saturn's moons keep the ring quasi-stable, which provides a 
critical threshold rate for collisional encounters due to the proximity of the co-orbiting 
ring particles. 

13.2.4 Magnetospheric Substorms and Auroras 

The size distribution of auroral areas have been measured with the UV Imager of the 
Polar spacecraft, which exhibits a powerlaw-like distribution with a slope of aA = 
1.21 ± 0.08 during active substorm time intervals, and aA — 1-16 ± 0.03 during 
quiescent time intervals (Fig. 13.6; Lui et al. 2000). The corresponding energy flux 
or power output P of auroral regions was derived to have a powerlaw slope of ap = 
1.05 ± 0.08 during active substorm time intervals, and ap = 1.00 ± 0.02 during 
quiescent time intervals (Lui et al. 2000). These powerlaw slopes are significantly 
flatter than predicted by our FD-SOC model, i.e., aA = 2.33 and ap = 1.67 for 
3D phenomena (S = 3, £'2 ~ 1-5). Although Lui et al. (2000) interpret auroras as 
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Figure 13.5: Top: Saturn's rings A, B, C, and the Cassini division, photographed by the Cassini 
spacecraft (credit: NASA, JPL, Space Science Institute). Bottom: Measurements of the particle 
size distribution functions for 8 ring regions with Voyager I radio occultation measurements 
(Ring C: CI. 35, CI. 51; Cassini division: CD2.01; Ring A: A2.12, A2.10, A2.14, A2.19, A2.24). 
The slopes of the fitted powerlaw functions in these 8 regions are: ol ~ 3.11, 3.05, 2.79, 2.74, 
2.70, 2.75, 2.93, 3.03. The range of particle sizes is L = 0.01 - 10 m (Zebker et al. 1985). 
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Fi gure 13.6! Left: Global image of the auroral oval observed by the Ultraviolet Imager (UVI) 
onboard the NASA satellite "Polar" on April 4, 1997 at 0519 UT, projected onto an Earth map 
(credit: NASA, Polar/UVI Team, George Parks). Right: Occurrence rate frequency distributions 
of auroral blobs as a function of the area (in units of square kilometers) during substorm-quiet 
time intervals, recorded with Polar LfVT during Jan 1-31, 1997 (Lui et al. 2000). 



a SOC phenomenon, the observed powerlaw slopes are far out of the range observed 
and predicted for other SOC phenomena. Thus, either the sampled distributions are 
incomplete, they underestimate the areas systematically for smaller events, or SOC 
models are not applicable for these events. 

Plasma flows in the magnetotail plasma with speeds v > 400 km s~^ were found to 
have a powerlaw distribution of durations T, with N{T) oc T^-59±o.07 (Angelopoulos 
et al. 1999), which is not too far off our theoretical prediction (with ar = 2.0), given 
the relatively small powerlaw range of only « 1.5 decades. Also the size distributions 
of the durations of AE index (a = 1.24; Takalo 1993; Takalo et al. 1999) and AU index 
(a — 1.3; Freeman et al. 2000; Chapman and Watkins 2001) were found to be much 
flatter than predicted by our SOC model. 

Electron bursts in the outer radiation belt (at 4 — 8 L-shell distances), which may be 
modulated by fluctuations of the solar wind, were found to have powerlaw distributions 
with slopes ofap = 1.5 — 2.1 and were interpreted as SOC phenomena (Crosby et 
al. 2005), which have slopes that are quite consistent with our SOC model (ap = 1.67). 
The solar wind is thought to be the source of these energetic electrons, although the 
solar wind velocity frequency distributions were found to exhibit significant deviations 
from simple powerlaws (Crosby et al. 2005). 
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Figure 13.7; Occurrence frequency distributions of hard X-ray peak count rates P [cts s ^] 
observed with HXRBS/SMM (1980 - 1989), BATSE (1991 - 2000), and RHESSI (2002 - 2010), 
with powerlaw fits. An average pre-flare background of 40 [cts s^^] was subtracted from the 
HXRBS count rates. Note that BATSE/CGRO has larger detector areas, and thus records higher 
count rates (Aschwanden 2011b). 

13.2.5 Solar Flares 

Solar flares are the best studied SOC phenomena in astrophysics. The impulsive energy 
release associated with solar flares, which can be observed in virtually all wavelengths, 
from gamma-rays, hard X-rays, soft X-rays, EUV, white-light, infrared, to radio wave- 
lengths, has been interpreted as a SOC phenomen from early on (Lu and Hamilton 
1991). Large datasets with n w 10'* — 10^ events sampled over up to eight orders of 
magnitude in energy provide the necessary statistics to determine accurate slopes of 
the observed powerlaw-like size distributions. However, major challenges exist still 
in the elimination of sampling biases in incomplete event sets, the understanding and 
modeling of powerlaw slopes in different wavelengths in terms of the underlying phys- 
ical scaling laws, and the automated determination of geometric parameters for large 
event datasets. A detailed account of observational results sorted into different wave- 
length regimes is given in Aschwanden (2012a; chapters 7 and 8). We summarize the 
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Figure 13.8: Occurrence frequency distributions of iiard X-ray total counts or fluence E [cts] 
observed with HXRBS/SMM (1980 - 1989), BATSE (1991 - 2000), and RHESSl (2002 - 2010), 
with powerlaw fits. An average pre-flare background of 40 cts s^^ multiplied with the flare 
duration was subtracted in the total counts of HXRBS (Aschwanden 2011b). 

results of powerlaw slopes observed in size distributions of various SOC parameters in 
Table 13.2, selecting mostly representative examples with large datasets from different 
instruments and wavelength regimes. 

Size distributions of peak fluxes P (Fig. 13.7), time-integrated fluxes or fluences 
E (Fig. 13.8), and flare durations T (Fig. 13.9), have been measured for energies > 
25 keV in hard X-ray wavelengths with instruments on the spacecraft ISEE-3, SMM, 
CGRO, and RHESSl, in soft X-ray wavelengths with Yohkoh and GOES, and in EUV 
with SOHO/EIT, TRACE, and AIA/SDO. Most of the observed powerlaw slopes were 
measured close to the theoretical predictions, i.e., ar = 2.0, ap — 1.67, and as = 
1.5 (Table 13.2), which is consistent with a dimensionality of 5* = 3, a mean fractal 
dimension of Ds ~ (1 + 5')/2 — 2.0, an energy-volume scaling index of 7 ss 1, and a 
diffusion power index of /3 « 1 (Eq. 13.17 and 13.19). 

The measurements in soft X-rays and hard X-rays are all made in broadband en- 
ergy and wavelength ranges, and thus are least biased regarding a complete sampling 
of all energy and temperature ranges. The probably most controversial measurements 
have been made for the smallest flares, also called nanoflares, which have typical tem- 
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Figure 13.9: Occurrence frequency distributions of hard X-ray flare durations T [s] observed 
with HXRBS/SMM (1980 - 1989), BATSE (1991 - 2000), and RHESSI (2002 - 2010) with 
powerlaw fits. The flare durations for RHESSI were estimated from the time difference between 
the start and peak time, because RHESSI flare durations were determined at a lower energy of 
12 keV (compared with 25 keV for HXRBS and BATSE), where thermal emission prolonges the 
nonthermal flare duration (Aschwanden 2011b). 



peratures of T « 1 — 2 MK and originate in small loops that barely stick out of the 
transition region. Since EUV measurements with SOHO/EIT, TRACE, and AlA/SDO 
are all made with narrowband temperature filters, the inferred thermal energies essen- 
tially scale with the flare area (Eq. 13.27) or volume (Eq 13.26), for which we predict 
powerlaw slopes of ath,A ~ ctA ~ 2.3 and ath,v ~ ctv ~ 2.0, which are significantly 
steeper than what is predicted for thermal energies sampled with broadband instru- 
ments, i.e., atfi w 1.5. This sampling bias has resulted into a controversy whether 
nanoflares dominate coronal heating, because a powerlaw slope steeper than the crit- 
ical value of 2 indicates that the energy integral diverges for the smallest events, as 
pointed out early on (Hudson 1991). Synthesizing measurements from narrowband 
EUV instruments with broadband soft X-ray instruments, as well as taking the fractal 
geometry of flare structures into account, however, could reconcile the size distribu- 
tion for nanoflares with that of large flares with a corrected value that is close to the 
theoretical prediction of as ~ 1-5 (Fig. 13.10; Aschwanden and Parnell 2002). 
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Figure 13.10: Synthesized frequency distributions from all three wavebands (TRACE 171 A, 
195 A, and Yohkoh/SXT AlMg) (grey histograms), along with the separate distributions from 
each waveband (in greyscales). Each of the distributions is fitted with a powerlaw, with the slope 
values and formal fit errors given in each panel. The four panels represent the four parameters 
of length L, area A, total emission measure EM (which is proportional to the peak flux P), and 
the thermal energy Eth (Aschwanden and Pamell 2002). 



Problematic are also the measurements of flare durations T for several reasons, such 
as the limited range of durations over which a powerlaw can be fitted, the ambiguity 
of separating overlapping long-duration flares, and the solar cycle dependence. While 
the event overlap problem is not severe during the solar minimum, where a slope close 
to the theoretically predicted value of aT = 2.0 is measured, the flare event pile-up 
bias becomes very severe during the solar maximum, producing powerlaw slopes of up 
to aT ^5 (Aschwanden and Freeland 2012). Also the powerlaw slope of hard X-ray 
peak counts ap appears to reveal a solar cycle dependence due to a similar effect (Bai 
1993, Biesecker 1994, Aschwanden 2011b). 

The least explored size distributions of solar flares are the length scale L and area A 
size distributions. Relatively small samples of flare areas have been measured for EUV 
nanoflares (Aschwanden et al. 2000, Aschwanden and Parnell 2002) and for the largest 
M and X-class flares (Aschwanden 2012b). Since the measurement of these parameters 
provides a direct test of the scale-free probability theorem (Eq. 13.1), without depend- 
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Figure 13.11: A solar flare event is observed on 2011 Mar 7, 19:43-20:58 UT, with AIA/SDO 
335 A, with GOES time profiles (top panel), the total EUV 335 A flux (second panel), the spatio- 
temporal evolution of the radius r{€) = ^J~MfypK of the time-integrated flare area A{t) for two 
thresholds, Fthresh = 100, 200 DN/s (third panel; histogrammed), fitted with the anomalous 
diffusion model (third panel; solid curve), and 5 snapshots of the baseline-subtracted flux (fourth 
row) and highpass-flltered flux (bottom row), with the threshold flux Fthresh = 100 DN/s shown 
as contour (Aschwanden 2012b). 



ing on any other physical parameter or model assumption, priority should be given 
to such measurements. Existing measurements have large error bars in the powerlaw 



13.2. OBSERVATIONS 



23 





20fl^y09 .' 


2^1^^^^ 


2^1^^ 


2011-|W5 


t 

2012^: 




it 




i0W.08-a9 


n 


2011-09-30 


2011-^^25^ 


2012-0; 


^06f 




2aii'r^-29 '' 




2011-10^22 






2'oim)2^ .• 


2011^6^ ■ 


201 1 :09l)7 


}|^^^^ 

201>v1JW^ 


2012W^23^"5^ 


2012-OC 








'20W19-1 




^^^^^^ 


\^ f- 

2012-01-27/ 


2012^3^^09^ 




2011-6,8.-03 \ 


201^^^ 






2^^^^(1: 


.■ .J- 
2012-03-10^- 


I:-.. - 


2011f0fe)4 




20fl-11-06 


2012-03-02 


2iff§§j 




2011-03-09 


t« 

2ai1-08H]8' ; 


1 

2011VR-24 




2011-11-09 




2012-03-15 



Figure 13.12: A selection of 48 solar GOES M and X-class flares observed with AIA/SDO at 
a wavelength of 335 A during 2010-2012. The contour shows the time-integrated flare area and 
the color-scale represents the intensity at 335 A. Note the complex spatial patterns that ressemble 
to the fractal avalanche patterns of cellular automaton simulations. 



slope due to the small number of analyzed events, but are largely consistent with the 
theoretical prediction of aL — 3.0 and aA — 2.33, calculated for 3D SOC avalanches 
(5" ^ 3) with a 2D (area S = 2) fractal dimension of D2 ~ {I + S)/2 ^ 1.5. 
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Besides the scale-free probability theorem, a second pillar of our FD-SOC model is 
the fractal-diffusive spatio-temporal relationship L{t) oc kT^/"^ (Eq. 13.8), which has 
been recently tested for a set of the 155 largest (GOES M and X-class) flare events (As- 
chwanden 2012b). An example of such a measurement of the spatio-temporal evolution 
is shown in Fig. 13.11, which exemplifies the spatial complexity of a large cluster of 
subsequent magnetic reconnection events that form a fractal volume filled with post- 
flare loops. Examples of another 48 large flares (GOES M and X-class) are shown 
in Fig. 13.12. Interestingly, the statistics of these 155 largest flares revealed a sub- 
diffusive regime (/? = 0.53 ± 0.27), while classical diffusion appears to be an upper 
limit. The diffusive characteristics measured in solar flares is consistent with the FD- 
SOC model, as it is was also found to be consistent with cellular automaton SOC 
simulations (Aschwanden 2012a). 

Radio bursts are produced in solar flares most frequently either by gyrosynchrotron 
emission of relativistic electrons that have been accelerated in magnetic reconnec- 
tion regions, or by electron beams that escape along magnetic field lines in upward- 
direction. Both types of radio bursts (microwave bursts and type III bursts) occur as 
a consequence of a plasma instability, and thus represent a highly nonlinear energy 
dissipation process that is typical for SOC processes. 

Somewhat ouf the predicted range are solar energetic particles, which have rather 
flat size distributions N{P) of the peak counts, typically in the range of ap « 1.1 — 1.4. 
However, since these are all high-energy particle events (> 10 MeV protons and > 3 
MeV electrons), we suspect that only the largest flares produce such high energies, and 
thus the sample is biased towards the largest flare events, which explains the flatter 
powerlaw slopes as a consequence of missing weaker events. 

13.2.6 Stellar Flares 

Impulsive flaring with rapid increases in the brightness in UV or EUV has been ob- 
served for a number of so-called flare stars, such as AD Leo, AB Dor, YZ Cmi, EK 
Dra, or e Eri. These types of stars include cool M dwarfs, brown dwarfs, A-type stars, 
giants, and binaries in the Hertzsprung-Russell diagram. Most of these stars are be- 
lieved to have hot soft X-ray emitting coronae, similar to our Sun (a G5 star), and thus 
magnetic reconnection processes are believed to operate in a similar way as on our Sun. 

However, what is different, is that the soft X-ray emission is several orders of mag- 
nitude stronger than from our Sun, if we put the stars into the same distance, and thus 
we expect an observational selection bias towards the largest possible flares. Another 
difference in flare statistics is that the meagre observational time allocation in the order 
of a few hours per star reveals only very few detectable flare events, in the order of 
5-15 events per observed star. A consequence of this small-number statistics is that 
we cannot determine a powerlaw slope from a log(N)-log(S) histogram, as we do for 
larger statistics (of at least 10^ up to 10''' events in solar flare data sets), but need to 
resort to rank-order plots, which correspond to the inverse distribution of cumulative 
occurrence frequency distributions. So, in principle, an inverse rank-order diagram can 
be plotted as shown in Fig. 13.13, which shows the logarithmic rank versus the flare en- 
ergy (or total counts) for each star, from which the powerlaw slope can be determined. 
However, if we deal with cumulative size distributions, we have also to be aware of the 
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Figure 13.13; Cumulative frequency distributions of flare energies (total counts) observed for 
12 cool (type F to M) stars with EUVE (Audard et al. 2000). The flare events are marked with 
diamonds, fitted with a powerlaw fit in the lower half (P; thick line), and fitted with a cumulative 
frequency distribution (C; ciu'ved function). 



drop-off that results at the upper end of the distribution due to the missing part in the 
powerlaw differential occurrence frequency distribution above the largest event. Thus, 
while a straight powerlaw function with slope a can be fitted to a differential frequency 
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distribution N{E), 

N{E) oc E-" , (13.32) 

the following function need to be fitted to the cumulative distribution Ncum{> E) 
(Aschwanden 2011a, section 7), 



J^^^- N{E')dE' jf'-E'-'-dE' {E^-'- - ElTa".) 

Ie::: ^iE')dE' j^—E'-<^dE' {e^-:: - E^-a-i) 

(13.33) 

The fit of this function to the cumulative distribution is shown in Fig. 13.13 for a set of 
12 flare stars, and the resulting values for the powerlaw slope of the inferred differ- 
ential occurrence frequency distributions (labeled with the letter C in Fig. 13.13). For 
comparison, we fit also a straight powerlaw with slope /3 = a — 1 to the lower half of 
the cumulative distribution, which is less effected by the upper cutoff (labeled with the 
letter P in Fig. 13.13). In Table 13.2 we summarize the means and standard deviations 
of the powerlaw slopes of flare energies observed on flare stars (see individual values 
in Table 7.7 of Aschwanden 201 la), which have been reported based on various other 
methods used by the authors, as = 2.17 ± 0.25. Fitting only the lower half of the 
distribution functions we find a significantly lower value of as = 1-85 ± 0.31, or by 
fitting Eq. (13.33) we find a similar range of as = 1-93 ± 0.35. This subtle differ- 
ence in the determination of the powerlaw slope is essential, because it discriminates 
whether the total energy radiated during stellar flares is dominated by the largest flares 
(if as < 2.0) or by nanoflares (if aE > 2.0; Hudson 1991). At this point it is not clear 
whether the difference in the powerlaw slopes obtained for stellar versus solar flares is 
due to a methodical problem of small-number statistics, or due to a sampling bias for 
super-large stellar flares, by solar standards. 



13.2.7 Pulsars 

Pulsars are fast-spinning neutron stars, which emit strictly periodic signals in radio 
wavelengths, as well as occasional giant pulses that represent glitches in the otherwise 
regular pulse amplitude and frequency. The glitches in pulse amplitude and frequency 
shifts correspond to large positive spin-ups of the neutron star, probably caused by 
sporadic unpinning of vortices that transfer momentum to the crust. Conservation of the 
angular momentum produces then an increase of the angular rotation rate. Thus, these 
giant pulses reveal highly nonlinear energy dissipation processes that can be considered 
as a SOC phenomenon and we expect a powerlaw function for their size distribution. 

Early measurements of the pulse height distribution of the Crab pulsar (NGC 0532 
or PSR B0531H-21) observed at 146 MHz were indeed found to have a powerlaw slope 
of ap ^ /3 + 1 — 3.5 over a range of 2.25 to 300 times the average pulse size, in a 
sample of 440 giant pulses (Argyle and Gower 1972). Similar values were measured 
by Lundgren et al. (1995) in a sample of 30,000 giant pulses, with ap « 3.06 — 3.36 
(Fig. 13.14, right). While the Crab pulsar is the youngest known pulsar (born in the year 
1054), PSR B1937H-21 is an older pulsar with a 20 times faster period (1.56 ms) than 
the Crab pulsar (33 ms). Cognard et al. (1996) measured a powerlaw distribution with 
a slope of ap « /3 -h 1 = 2.8 ± 0.1 from a sample of 60 giant pulses for this pulsar. 



13.2. OBSERVATIONS 



27 




Log f (Jy) 



Figure 13.14: Left: Crab nebula, which harbors the pulsar in the center (photographed by 
Hubble Space Telescope, NASA). Right: Frequency distribution of giant-pulse tlux densities 
measured from the Crab pulsar, observed during 15-27 May 1991 with the Green Bank 43-m 
telescope at 1330, 800, and 812.5 MHz. The tail can be represented by a powerlaw distribution 
Nf oc with a = 3.46 ± 0.04 for fluxes F > 200 Jy (Lundgren et al. 1995). 

Statistics of nine other pulsars revealed powerlaw slopes in a much larger range of 
aE = —0.13, ...,2.4 for the size distribution of pulse glitches (Melatos et al. 2008), but 
those measurements were obtained from much smaller samples of 6-30 giant pulses, 
and thus represent small-number statistics. 

The typical value of ap « 3.0 found in two pulsars deviates significantly from 
the prediction (ap = 1.67) of our FD-SOC model, and thus requires either a different 
model or more statistics from other pulsar cases. Preliminary values of small samples 
from other pulsars indicate a wide range (as = —0.13, 2.4, Melatos et al. 2008) 
that do not point towards a particular value that can be explained with a single model. 

13.2.8 Soft Gamma-Ray Repeaters 

Gamma-ray bursts were observed from a variety of astrophysical objects, such as neu- 
tron stars or black holes, but usually only one burst has been observed from each object. 
An exception is a class of objects that show repetitive emission at low-energy gamma- 
rays (> 25 keV), which were termed soft gamma-ray repeaters (GRS). Observations 
with the Compton Gamma Ray Observatory (CGRO) revealed four such SGR sources 
up to 1999, three in our galaxy and one in the Magellanic Cloud). At least three of these 
SGR objects were associated with slowly rotating, extremely magnetized neutron stars 
that are located in supernova remnants (Kouveliotou et al. 1998, 1999). It is believed 
that these soft gamma-ray bursts occur from neutron star crust fractures driven by the 
stress of an evolving, ultrastrong magnetic field in the order of i? > 10^* G. 

Occurrence frequency distributions of the fluence of soft gamma-ray repeaters were 
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Figure 13.15: Left: The Soft Gamma-Ray Repeater SGR 1900+14 in the constellation of 
Aquila photographed with the Spitzer Space Telescope in infrared (courtesy: NASA). Right: 
Differential frequency distributions of the fluences of soft gamma-ray repeater SGR 1900+14 
observed with CGRO, RXTE, and ICE (Gogus et al. 1999). 



obtained from four SGR sources: a database of 837 gamma-ray bursts from SGR 
1900+14 during the 1998-1999 active phase showed a powerlaw slope of as = 1-66 
over 4 orders of magnitude (Fig. 13.15; Gogus et al. 1999); and a combined database 
from SGR 1806-20, using 290 events detected with the Rossi X-Ray Timing Explorer, 
111 events detected with CGRO/BATSE, and 134 events detected with the Interna- 
tional Cometary Explorer (ICE), showing power laws with slopes of a£;=1.43, 1.76, 
and 1.67 (Gogus et al. 2000). These measurements agree remarkably well with the 
frequency distributions predicted by the FD-SOC model (ue — 1.50) as well as with 
those observed during solar flares, which were also observed at the same hard X-ray 
energies of > 25 keV. However, the physical energy dissipation mechanism may be 
quite different in a solar-like star and a highly magnetized neutron star, given the huge 
difference in magnetic field strengths [B w 10^ G for solar flares versus B « 10^** 
G in a magnetar), although magnetic reconnection processes could be involved in both 
cases. Nevertheless, soft gamma-ray repeaters have been interpreted as a SOC sys- 
tem (Gogus et al. 1999), in terms of a neutron star crustquake model (Thompson and 
Duncan 1996), in analogy to the SOC interpretation of earthquakes. 

13.2.9 Black Hole Objects 

Cygnus X-1, a galactic X-ray source in the constellation Cygnus, is the first X-ray 
source that has been widely accepted to be a black-hole candidate. The mass of Cygnus 
X-1 is estimated to be about 14.8 solar masses and it has been inferred that the object 
with (an event horizon at) a radius of 26 km is far too compact to be a normal star. 
Cygnus X-1 is a high-mass X-ray binary star system, which draws mass from a blue 
supergiant variable star (HDE 226868) in an orbit of 0.2 AU around the black hole. 
The stellar wind of this blue companion star swirls mass onto an accretion disk around 



13.2. OBSERVATIONS 



29 




Figure 13.16: Left: This artist concept of Cygnus X-1 shows the black hole drawing material 
from a companion star (right) into a hot, swirling accretion disk that surrounds the (invisible) 
black hole (Courtesy of Chandra X-Ray Observatory, NASA). Right: Observed frequency dis- 
tribution of the peak intensities of pulses in the light curve of the black-hole object Cygnus X-1, 
exhibiting a powerlaw slope of ap ~ 7.1 (Negoro et al. 1995; Mineshige and Negoro 1999). 

the black hole (Fig. 13.16). The X-ray time profile from Cygnus X-1 reveals time 
variability down to 1 ms, which is attributed to X-ray pulses from matter infalling 
toward the black hole and the resulting turbulence in the accretion disk. 

Observations of the X-ray light curve of Cygnus X-1 with Ginga exhibit a complex 
power spectrum that entails at least 3 piece-wise powerlaw sections, which have been 
interpreted as a superposition of multiple 1/f-noise spectra (Takeuchi et al. 1995). The 
occurrence frequency distribution of peak intensities shows a powerlaw-like function 
with a steep slope of ap « 7.1 (Fig. 13.16 right; Negoro et al. 1995; Mineshige and 
Negoro 1999). All these properties have been modeled with a sophisticated cellular 
automaton model in the framework of the SOC concept. Infalling mass lumps in the 
accretion disk are thought to trigger turbulent instabilities in the neighborhood of an in- 
fall side, which propagate avalanche-like and produce hard X-rays either by collisional 
bremsstrahlung or some other magnetically driven instability (e.g., a magnetic recon- 
nection process). The cellular automaton simulations (Takeuchi et al. 1995, Mineshige 
and Negoro 1999) were able to reproduce steep powerlaw slopes of the peak fluxes in 
the range of ap w 5.6 — 11.5, depending on the effect of enhanced mass transfer by 
gradual diffusion in addition to the avalanche-like shots, and this way could reproduce 
the observations. 

We note that the observed steep powerlaw slopes of peak fluxes (ap « 7.1) exceed 
the predictions of the FD-SOC model (ap = 1.67) by far. Such a steep slope, ap = 
1 + 2/(87) ~ 7, can be produced by an extremely weak dependence of the X-ray 
flux F on the avalanche volume V, i.e., F oc V with 7 « 1/9, which is different 
from the flux-volume scaling law of optically-thin soft X-ray or EUV emission (7 w 
1) generally observed in astrophysical sources, and thus may indicate a nonthermal 
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Figure 13.17: Left: The anatomy of a blazar is shown in a top view and side view, having 
the rotation axis with the relativistic jet oriented toward Earth or orthogonal to the line-of-sight 
(Courtesy of The Encylopedia of Science). Right: Frequency distribution of peak fluxes of flar- 
ing events in blazar GC 0109+224, including fluxes above a 3(7-threshold, fitted with a powerlaw 
function N{P) oc p-i-s^ (Ciprini et al. 2003). 

emission mechanism. Alternatively, the very steep powerlaw slope of ap could by 
part of an exponential frequency distribution near the upper cutoff, which could only 
be proven by sampling peak fluxes with higher sensitivity. Simultaneous modeling of 
the observed occurrence frequency distributions of time scales, peak fluxes, and power 
spectra in terms of a SOC model may reveal the underlying physical scaling law of the 
emission mechanism in black-hole accretion disks. 

13.2.10 Blazars 

Blazars (blazing quasi-stellar objects) are a special subgroup of quasars. This group 
includes BL Lacertae objects, high polarization quasars, and optically violent variables. 
They are believed to be active galactic nuclei whose jets are aligned within < 10° of 
our Une-of-sight (Fig. 13.17 left). Because we are observing "down" the jet direction 
we observe a large degree of variability and apparent superluminal speeds from the 
jet-aligned emission. The structure of blazars, like all active galactic nuclie, is thought 
to be powered by material falling onto the supermassive black hole at the center of 
the host galaxy, which emits highly intermittent gyrosynchrotron emission (in radio 
wavelengths), inverse Compton emission (in X-rays and gamma-rays), and free-free 
bremsstrhalung emission (in soft X-rays), modulated by the variable rate of matter 
infalling into the accretion disk of the central black hole. 
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Figure 13.18: Left: The diagram shows the cascade shower of a primary cosmic ray particle, 
which produces in the Earth's atmosphere a shower of secondary particles, which can be detected 
by an unterground detector, such as Cosmo ALEPH (courtesy of Cosmo ALEPH Team). Right: 
Cosmic ray spectrum in the energy range of iJ = 10^ — 10^^ eV, covering over 12 orders of 
magnitude. There is a "knee" in the spectrum around E ~ 10^'' eV, which separates cosmic rays 
originating within our galaxy (at lower energies) and those from outside the galaxy (at higher 
energies). (Credit: Simon Swordy, University of Chicago). 

The optical variability of blazar GX 0109+224 was monitored and the light curves 
were found to exhibit flickering and shot noise, with a power spectrum P(i^) oc i/^^ 
with power index of p = 1.57 — 2.05 (Ciprini et al. 2003). The occurrence frequency 
distribution of peak fluxes of flare events was found to have a powerlaw slope of ap — 
1.55 (Fig. 13.17, right; Ciprini et al. 2003), which is close to the prediction of the FD- 
SOC model (ap — 1.67). Thus, the highly- variable blazar emission was interpreted 
in terms of SOC models. The fact that the peak size distribution of radio emission 
observed in the blazar agrees with the prediction of the FD-SOC model is consistent 
with a near-proportional radio flux-volume scaling, i.e., F oc y with 7 ~ 1, which is 
generally the case for gyro-synchrotron emission. This is different from the flux scaling 
of emission observed from the black hole Cygnus X-1. Thus, SOC statistics allows us 
to discriminate between different physical emission mechanisms in black holes and 
blazars. 

13.2.11 Cosmic Rays 

Cosmic rays are high-energy particles that have been accelerated during a long journey 
through a large part of our universe, inside our Milky way as well as from outside of 
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our galexy. Their boosting to the highest energies of up to < 10^^ eV can only oc- 
cur by multi-step acceleration processes throughout the universe over long time spans. 
The energy spectrum of cosmic rays, as it can be measured from showers of secondary 
particles produced in the Earth atmosphere (Fig. 13.18 left), exhibits a powerlaw-like 
energy spectrum extending from 10^ eV (=1 GeV) over 12 orders of magnitude up 
to < 10^^ erg (Fig. 13.18 right). The average slope is as ~ 2.7. However, a more 
accurate model is a broken powerlaw with a "knee" in the spectrum around E w 10^^ 
eV. The widely accepted interpretation of this knee is that it separates the origin of cos- 
mic rays from inside and outside of our galaxy. The powerlaw slope above the "knee" 
steepens to aE ~ 3.0 — 3.3. Interestingly, our FD-SOC model applied to the kinetic 
energy gain in a coherent direct current (DC) electric field implies a proportionality of 
El (X L (Eq. 13.28) for sub-relativistic energies, and thus predicts an energy distribu- 
tion of N{El) cx EJ^^, similar to the cosmic ray spectrum. Of course, cosmic rays 
are highly relativistic and are likely to be produced by many (n) acceleration phases. 
However, even if each acceleration phase is local and has a relatively small length of 
/ ^ L, the energy gain of the particle would add up linearly with increasing travel time 
and travel distance J2 ' ^^'^ could still fulfill the proportionality. 



and end up with an energy spectrum of N{El) oc L~^. 

So, can we understand the acceleration of a cosmic ray as a SOC process? The 
lattice grid would cover a large fraction of the extragalactic space of our universe, 
the ensemble of cosmic-ray particles would represent an avalanche that nonlinearly 
dissipates energy from local acceleration processes (such as by the first-order Fermi 
acceleration process between intergalactic or interstellar magnetic clouds), which are 
self-organizing in the sense that accelerating fields (of magnetic clouds) are constantly 
restored by the galactic and interstellar dynamics. 

Our fractal-diffusive spatio-temporal relationship (r{t) oc t^^"^, Eq. 13.8) further 
predicts a random- walk through the universe and an age T that scales as L oc T^/^ with 
the straight travel distance L from the point of origin of the cosmic ray particle. Using 
the "knee" in the cosmic-ray energy spectrum (at Egai ~ 10^^ eV) as a calibration 
for the distance of the Earth to the center of our galaxy (Lgai » 50 fight years ~ 
5 X 10^^ cm), we can estimate the straight length scale over which the cosmic ray 
particle travelled by random walk 



which yields L^ax « 5 x 10^ erg, which corresponds to about 10% of the size of 
our universe r„„i « 4 x 10^^. Since the acceleration efficiency is different in galactic 
and extragalactic space, the diffusion coefficient of the random walk is also different 
and therefore we expect a different powerlaw slope in the energy spectrum produced in 
these two regimes. 



n 



n 




(13.34) 




(13.35) 
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13.3 Conclusions 

In this chapter we generalized the fractal-diffusive self-organized criticality (FD-SOC) 
model in terms of four fundamental parameters: (i) the Euclidean dimension S, (ii) the 
fractal dimension Ds of the spatial SOC avalanche structure, (iii) the diffusion index /3 
that includes both sub-diffusion and super-diffusion, and (iv) the energy-volume scal- 
ing law with powerlaw index 7. This model predicts powerlaw functions for the occur- 
rence frequency distributions of the SOC model, and moreover predicts their powerlaw 
slope as a function of the four fundamental parameters. For a Euclidean dimension of 
S = 3,a mean fractal dimension of Ds f« (1 + S')/2 = 2, classical diffusion (/3 = 1), 
and linear flux-volume scaling (7 = 1), our generalized FD-SOC model predicts then 
the following powerlaw slopes: a/, = 3 for length scales, ua = 2.333 for areas, 
aT = 2 for durations, a = 1.667 for peak fluxes, and as = 1-5 for fluences or total 
energies of the SOC avalanches. 

Comparing these theoretical predictions with the observed powerlaws of size distri- 
butions in astrophysical systems (summarized in Table 13.2) we find acceptable agree- 
ment for the cases of lunar craters, asteroid belts, Saturn rings, outer radiation belt 
electron bursts, solar flares, soft garmna-ray repeaters, and blazars, if we apply the lin- 
ear flux- volume scaling. Discrepancies are found for magnetospheric substorms, steUar 
flares, pulsar glitches, black holes, and cosmic rays, which apparently require a nonlin- 
ear flux- volume scaUng. Pulsar glitches and cosmic rays can indeed be modeled by as- 
suming a linear energy-length scaling, which leads to energy spectra of N{E) ^ E~^. 
Black-hole pulses have very steep size spectra, which indicates a quenching or satura- 
tion process that prevents a large variation of pulse amplitudes. Magnetospheric sub- 
storms and solar energetic particles have the flattest size distributions, which possibly 
can be explained by a selection effect with a bias for the largest events. In conclusion, 
the generalized FD-SOC model can explain a large number of astrophysical observa- 
tions and can discriminate between different scaUng laws of astrophysical observables. 
We envision that more refined scaUng laws between astrophysical observables will be 
developed that are consistent with the observed size distributions, and this way will 
provide the ultimate predictive power for SOC models. 
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Table 13.2: Summary of theoretically predicted and observed powerlaw indices of size distri- 
butions in astrophysical systems. 





Length 


Area 


Duration 

aT 


Peak flux 

ap 


Fluence 


FD-SOC Theory 


3.0 


2.33 


2.0 


1.67 


1.50 


Lunar craters: 

Mare TranquilUtatis ^) 

Meteorites and debris ^) 


3.0 
2.75 










Asteroid belt: 
Spacewatch Surveys^) 
Sloan Survey'^) 
Subaru Survey') 


2.8 

2.3-4.0 
2.3 










Saturn ring: 
Voyager 1®) 


2.74-3.11 










Magnetosphere: 
Active substorms^) 
Quiet substorms^) 
Substorm flow bursts*) 
AE index bursts®) 
AU index bursts^") 
Outer radiation belt^^) 




1.21 ±0.08 
1.16 ±0.03 


1.59 ±0.07 

1.24 

1.3 


1.05 ±0.08 
1.00 ±0.02 

1.5-2.1 




Solar Flares: 
lSEE-3, HXR12 
HXRBS/SMM, HXR13 
BATSE/CGRO, HXR^^ 
RHESSl, HXR15 
Yohkoh, SXR16 
GOES, SXRi^ 
SOHO/EIT, EUV^^ 
TRACE, EUyi® 
AIA/SDO, 335 A, EUV^^ 
Microwave bursts^"*^ 
Type 111 bursts22 
Solar energetic particles^^ 


1.96-2.41 

2.50-2.75 
1.96 


1.77-1.94 

2.3- 2.6 

2.4- 2.6 


1.88-2.73 
2.17 ±0.05 
2.20-2.42 
1.8-2.2 

2.0-5.0 
1.4-2.0 

2.17 


1.75-1.86 

1.73 ±0.01 

1.67-1.69 

1.58-1.77 

1.64-1.89 

1.86-1.98 

1.52-2.35 

1.34 

1.2-2.5 

1.26-1.91 

1.10-2.42 


1.51-1.62 

1.53 ±0.02 

1.56-1.58 

1.65-1.77 

1.4-1.6 

1.88 

1.41-2.06 
1.27-1.32 


Stellar Flares: 
Flare stars (reported)^* 
Flare stars (powerlaw fit)^^ 
Flare stars (cumulative fit)^^ 










2.17 ±0.25 
1.85 ±0.31 
1.93 ±0.35 


Astrophvsical Objects: 
Crab pulsar^^ 
PSRB1937+2126 

Soft Gamma-Ray repeaters^*^ 
Cygnus X-1 black hole^* 
Blazar GC 0109H-22429 
Cosmic rays^° 








3.06-3.50 
2.8 ±0.1 

7.1 
1.55 


1.43-1.76 
2.7-3.3 
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